Kondo effect in coupled quantum dots with RKKY interaction: 
Finite temperature and magnetic field effects 
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We study transport through two quantum dots coupled by an RKKY interaction as a function of 
temperature and magnetic field. By applying the Numerical Renormalization Group (NRG) method 
we obtain the transmission and the linear conductance. At zero temperature and magnetic field, 
we observe a quantum phase transition between the Kondo screened state and a local spin singlet 
as the RKKY interaction is tuned. Above the critical RKKY coupling the Kondo peak is split. 
However, we find that both finite temperature and magnetic field restore the Kondo resonance. Our 
results agree well with recent transport experiments on gold grain quantum dots in the presence of 
magnetic impurities 

PACS numbers: 73.63.Kv, 72.15.Qm, 71.27.+a, 73.23.Hk 



Introduction. In recent years the Kondo effect [l| in 
semiconductor quantum dots has gained significant in- 
terest both theoretically and experimentally |M 
Electronic transport in quantum dots is strongly influ- 
enced by the Coulomb blockade (^J due to their small 
size. In these systems, a single unpaired spin can inter- 
act with conduction electrons, leading to screening of the 
spin and an enhanced conductance at low bias and low 
temperatures. More recently, due to rapid progress in 
spintronics and quantum information, it becomes desir- 
able to gain more tunable spin control in double quantum 
dots where an effective spin-spin interaction known as 
the Ruderman-Kittel-Kasuya-Yoshida (RKKY)P cou- 
pling is generated between the two dots via conduction 
electrons in the leads. The RKKY coupling competes 
with the Kondo effect in these systems, leading to a 
quantum phase transition between the Kondo screened 
phase at weak RKKY coupling and a local spin-singlet 
state at strong coupling. Part of this rich physics has 
been studied previously in the two-impurity Kondo 
and Anderson 8] problems. 

Recently, indications for a competition between Kondo 
screening and the RKKY interaction have been observed 
experimentally 0. This experiment stimulated theo- 
retical and experimental efforts on the above-mentioned 
quantum phase transition in coupled dots. Very re- 
cently, new measurements have been reported to observe 
the restoring of the Kondo resonance in a gold quan- 
tum dot with magnetic impurities at finite temperatures 
and magnetic fields where an effective RKKY cou pling 
is generated between the dot and the impurities |lfll |. 
Though similar systems have been studied theoretically 
via several approaches [ill IT3. Il3j | , little is known about 
the transport properties at finite temperatures and mag- 
netic fields. Motivated by these recent experiments, here 
we systematically study transport properties of double 
quantum dots coupled by an RKKY interaction at fi- 
nite temperatures and magnetic fields via the Numeri- 



cal Renormalization Group (NRG) method [lj|, a non- 
perturbative approach to quantum impurity systems. In 
contrast to other techniques, this method does not rely 
on any assumptions regarding the ground state or the 
leading divergent couplings, which is crucial in our anal- 
ysis. 

The Model. We consider a two-impurity Anderson 
model as shown in Fig. ^ describing two quantum dots 
which are seperately coupled to two-channel leads and 
subject to an antiferromagnetic RKKY spin-spin inter- 
action and a local magnetic field. This setup is general 
enough to describe both experiments mentioned above. 
The Hamiltonian of the system is given by: 
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Here, a — L/R denotes the left/right lead, i — 1,2 de- 
notes the two dots, Ni = Y^^ dl„d- is the number of elec- 
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trons occupying dot i and = {^I'^Mlaa' A(7^<y<'''^ia' 
are the spins of the two levels. Each dot is subject to a 
charging energy U and J is the RKKY exchange coupling 
between the two dots which is assumed to be antiferro- 
magnetic (J > 0). The magnetic field B induces a Zee- 
man splitting, where we have set g = 1. Here, we consider 
the model with particle- hole symmetry, i.e. e^i = — -^i 
and single occupation for each dot, i.e. M ~\. The en- 
ergies e^i of the two dots and their precise position in the 
Coulomb blockade valley can be tuned experimentally 
by an external magnetic field and the gate voltage. Here, 
we neglect the energy dependence of the tunneling ma- 
trix elements Via ■ Without loss of generality, we assume 



2 



they have left-right symmetry, i.e. Vil = Vm = Vi, 
V2L = V2R = V2, but Vi ^ V2 in general. The intrinsic 
line width of the dot levels due to tunnel coupling to the 
leads is r = + Tfl with r^/fl = 2'k\V\'^Nl/i^ where 
N[^/f( is the density of states in the leads. 



with the Fermi function f{uj) and the transmission coef- 
ficient (IJ] 
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FIG. 1: Two quantum dots (denoted by the two circles) with 
two-channel leads (chemical potentials fiL and hr). Vi and 
V2 denote the tunneling matrix element between the dots and 
the leads. There exists an RKKY coupling J between the two 
dots. 

For J = 0, and Vir = 0, three triplet configurations 
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^^|0), |1,0) - (1/V2)(4^4^ +4^4^10), 

= 4i4il0) and the singlet |0,0) = dlt^xlO) 
are degenerate. Finite tunneling Vir into the leads to 
independent spin- 1/2 Kondo screening in each of the 
two dots. The corresponding Kondo temperatures 
and are generally different, and given by T^^'' cc 
D exp{—l/2pcJi(2)) where D is the bandwidth, pc the 
density of states of the leads, and Ji(2) = 4F^2)/^ 
effective Kondo coupling 0. 

The above degeneracy is lifted at finite RKKY coupling 
J: three triplet states are shifted to energy Et — J/4 and 
the singlet state to energy Eg = — 3J/4. There exists 
competition between Kondo screening and a local spin- 
singlet ground state: the former is expected to be the 
ground state for large J, the latter for small J. A quan- 
tum phase transition at zero temperature between these 
two phases occurs as the RKKY coupling is tuned 0, ■ 
Note that a related singlet-triplet transition in two-level 
quantum dots was studied both for two-mode and 
for single-mode leads 0, 0, 0|. In the latter case, a 
Kosterlitz-Thouless transition from a local singlet state 
to a single-channel 5=1 underscreened Kondo model 
was observed [lij . 

An interesting aspect of the Kondo-to-singlet transi- 
tion in our setup is the tunability between these two 
phases since one can get good control over the various 
parameters in experiments. In particular, as observed 
in the experiment [loj . the Kondo resonance is restored 
at finite temperatures and magnetic fields close to the 
singlet-triplet transition. The goal of our work is to de- 
scribe this behavior theoretically. 

Transport properties. We are interested in calculating 
electronic transport through the dot close to the tran- 
sition. To this end, we use the generalized Landauer 
formula 
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Here we have introduced the retarded dot Green's func- 
tions Gii'ait) = -iO{t){{d^^{t),dl,J). In the following 
we focus on the low bias regime, where T(w) can be eval- 
uated in equilibrium, using the NRG. For a detailed de- 
scription of this technique see Ref . [l3| . Using the current 
formula (jSJ, we determine the behavior of the linear con- 
ductance G{T) — ^\v=o at finite temperatures. Note 
that the equilibrium transmission T(uj) also yields a good 
approximation to the differential conductance dl/dV at 
finite bias measured in experiments. 
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FIG. 2: (a) Transmission coefficient at zero temperature for 
different RKKY interactions. The unit of energy is the half 
bandwidth D = 1 of the conduction electrons. Here, U = 1, 
Edl = <^d2 = —0.5 , FlL = Flij = 0.05, = r2fl = 0.1, 



A = 



0.002, 



/(c.-Mfl))TH (2) 



0.00002 (for J = 0). The critical 
RKKY coupling is Jc ~ 0.00042. (b) Linear conductance 
G(T) for different RKKY couplings. Parameters are the same 
as in (a). 



Quantum phase transition at zero temperature for zero 
field. In Fig. |2Ia) we plot the transmission coefficient 
T{w) as a function of RKKY coupling J for zero temper- 
ature and zero field. It shows a quantum phase transition 
between the Kondo screened phase (J < Jc) and the lo- 
cal spin-singlet phase (for J > Jc), where < Jc <T^. 
The value of Jc is consistent with previous results by 
Sakai et al. 8]. For J < Jc, a Kondo peak at w = is ob- 
served, where the transmission T(w — s- 0) tends to reach 
its unitary limit of 2 (each dot acts as an unitary Kondo 
channel). Due to systematic numerical errors in the NRG 
calculation, this limit is underestimated by about 10%. 
As J is increased, the Kondo peak becomes narrower, 
indicating a vanishing low-energy scale as the system ap- 
proaches the critical point. As J ^ Jc, the transmission 
reaches the unstable fixed point value T{u! ^ 0) = 1. 
For J > Jc, the two dots form a local spin-singlet state 
where the Kondo peak is split and T{u! ^ 0) = 0. The 
splitting of the Kondo peak increases as J is increased. 

Transport at finite temperatures. We now determine 
the behaviour of the linear conductance G{T)/Gq at fi- 
nite temperatures where Gq = 2e^ /h is the conductance 
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unit. Results are shown in Fig. |2Ib). The Kondo screened 
and local spin-singlet phases are characterized by stable 
fixed points with G{T — > 0) « 2 or 0, respectively. In the 
Kondo regime, upon lowering the temperature, the con- 
ductance rises up to the unitary limit with two step-like 
structures indicating the two crossover Kondo tempera- 
tures and r|. 

In the local spin-singlet regime we find a nonmonotonic 
behaviour of the conductance when T is lowered: After 
an initial rise due to the Kondo effect, G{T) decreases to 
zero as the temperature is lowered. The broad peak of 
G{T) at finite temperatures is in fact a signature of the 
reappearance of the Kondo resonance which is also seen 
in the experiment |inj |. 



T(co) 




FIG. 3: Transmission coefficient r(w) for a fixed RKKY cou- 
pling J — 0.002 at finite temperatures. Here, Til = Tir — 
0.04, T2L ~ = 0.08. The remaining parameters are the 
same as in Fig. 2. 

To gain more insight into the finite temperature behav- 
ior, we present the plot of transmission T(u)) at different 
temperatures (see Fig. O when RKKY is strong enough 
(J > Jc) so that the Kondo peak is split at zero tempera- 
ture. We find that as temperature increases, r(w) inside 
the dip first increases due to thermal broadening of the 
split peaks around w ±J until the Kondo resonance 
reappears. Then the peak height decreases again as T 
is further increased, similar to the Kondo effect at finite 
temperatures without RKKY interaction. The transmis- 
sion T{uj = 0) reaches its maximum value at a tempera- 
ture Tmax oc J — Jc- Note that the NRG data for T{ui) at 
Lo <T are estimated by extrapolating T{u}) from higher 
frequencies. The qualitative behavior of the transmission 
- reappearance of the Kondo peak at finite T - agrees well 
with recent experiments jl(Tj |. 

Effect of a magnetic field. In the presence of a mag- 
netic field a singlet-triplet crossover occurs. Here we dis- 
cuss the effect of a Zeeman splitting at zero temperature 
in the presence of a large RKKY coupling J such that 
the Kondo peak is split in the absence of magnetic fields 
(see Fig. 0J. For finite B < J, wc find an increase of 
the transmission T{u!) a,t ui — 0. When B is compara- 
ble to the value of RKKY interaction, B ^ J, we ob- 
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FIG. 4: Transmission coefficient T{w) of the coupled dots at 
different magnetic fields. Here, Fil = Tir = = ^2r = 
0.06, J = 0.002. The remaining parameters are the same as 
in Fig. 121 Inset (a): the transmission ^(a;) for B = 0.0016, 
which clearly shows four peaks in pairs around ±(B ± J). 
Inset (b): The energy levels of the triplet and singlet states in 
the presence of a finite RKKY coupling J and a magnetic field 
B. The three triplet states are split into = 1/4 J — B, 
Ex,o = 1/4J, and Ex,-i = 1/4J -|- B. The singlet |0, 0) and 
one of the triplet states jl, 1) become degenerate at B = J. 



serve the reappearance of the Kondo peak where T(uj) 
reaches the unitary limit of 1, corresponding to a single 
channel 5=1/2 Kondo effect. When the magnetic field 
increases further, the Kondo peak splits again. What 
happens is that due to the Zeeman splitting, one compo- 
nent of the triplet (|1, 1)) is "pulled down" and eventu- 
ally becomes degenerate with the singlet |0,0) (see inset 
(b) of Fig. 21). At this point a Kondo effect - analogous 
to S = 1/2 Kondo - arises between these two states, 
which has been discussed previously for a 2-level quan- 
tum dot within a perturbative scaling approach |20l | . Our 
nonperturbative NRG results confirm this scenario: The 
transmission T(w = 0) indeed reaches the unitary limit 
of the S* = 1/2 Kondo effect. Note that due to tunnel- 
ing into the leads, degeneracy of the singlet and triplet 
state occurs at a renormalized value B ^ J . Close to the 
degeneracy point, a four-peak structure emerges in the 
transmission (see inset (a) of Fig. ^J, corresponding to 
the splitting between the singlet and the lowest (second- 
lowest) triplet state. Due to NRG broadening effects, the 
third triplet state is not visible. 

In Fig. 13 we present results for the singlet-triplet 
crossover at smaller RKKY coupling, for different mag- 
netic fields and finite temperature. Compared to the 
sharp transition in Fig. ^ where J is much larger than 
Tfe, the crossover here is smoother and closer to the ex- 
perimental observation in Ref. [l3- Note that due to ther- 
mal broadening at a finite temperature, T(u! = 0) > 
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even in the absence of a magnetic field. For antiferromag- 
netic RKKY coupling (see Fig.[Sfa)), T{uj) shows the dip- 
peak-dip structure with increasing magnetic field. For 
ferromagnetic RKKY (see Fig.lSfb)), the Kondo peak at 
B = 0, which is due to complete screening of the triplet, 
splits monotonically with increasing field due to Zeeman 
splitting of the triplet levels as shown in Ref. ^3- Both 
results are qualitatively in good agreement with experi- 
ment 113 (for new measurements on the magnetic field 
dependence see H^l)- 
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FIG. 5: (a) Transmission coefficient T{uj) of the coupled 
quantum dots for different magnetic fields. Here, U — 1, 
Edi = ed2 = -0.5 , FiL = Fiji = = r2_R = 0.1, A = 4, 
J = 0.007, Tk ~ 0.0025 (for J = 0), Jc^ 0.005, T = 0.00001. 
The values of B are in steps of 0.001, and the traces of T(ijj) 
are shifted in steps of 4005. (b) Transmission vs. frequency 
for ferromagnetic RKKY coupling J — —0.005 and different 
magnetic fields. The remaining parameters are the same as 
in (a). 

Conclusions. We have studied transport properties 
of a double quantum dot system with RKKY interac- 
tion. Using the Numerical Renormalization Group we 
have calculated the transmission at finite frequency as 
a function of temperature and magnetic field. A quan- 
tum phase transition between the Kondo screened phase 
and a local spin-singlet state is observed. Moreover we 
have shown that both finite temperature and a magnetic 
field can restore the Kondo resonance in the presence of 
an RKKY coupling. This crossover back into a Kondo 
screened state is in good agreement with recent measure- 
ments |lOj . For the differential conductance in a finite 
magnetic field we predict a multiple peak structure which 
is yet to be observed in future experiments. 
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